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In astro-physical environments with intense neutrino fluxes, neutrino-neutrino forward scattering 
contributes both diagonal and off-diagonal potentials to the flavor-basis Hamiltonian that governs 
neutrino flavor evolution. We examine a special case where adiabatic flavor evolution can produce 
an off-diagonal potential from neutrino-neutrino forward scattering that dominates over both the 
corresponding diagonal term and the potential from neutrino-matter forward scattering. In this 
case, we find a solution that, unlike the ordinary Mikeyhev-Smirnov-Wolfenstein scenario, has both 
neutrinos and antineutrinos maximally mixed in medium over appreciable ranges of neutrino and 
antineutrino energy. Employing the measured solar and atmospheric neutrino mass-squared differ- 
ences, we identify the conditions on neutrino fluxes that are required for this solution to exist deep 
in the supernova environment, where it could affect the neutrino signal, heavy-element nucleosyn- 
thesis, and even the revival of the supernova shock. We speculate on how this solution might or 
might not be attained in realistic supernova evolution. Though this solution is ephemeral in time 
and/or space in supernovae, it may signal the onset of subsequent appreciable flavor mixing for both 
neutrinos and antineutrinos. A similar solution may also exist in an early universe with significant 
net neutrino-lepton numbers. 

PACS numbers: 97.60.Gb,14.60.Pq 



I. INTRODUCTION 

In this paper we study the problem of coherent nonlin- 
ear flavor evolution of active neutrinos in environments 
where neutrino-neutrino forward scattering makes a sig- 
nificant contribution to the effective neutrino mass in 
medium. In particular, we examine a special case where 
the off-diagonal potential from neutrino-neutrino forward 
scattering becomes the dominant term in the flavor-basis 
Hamiltonian that governs neutrino flavor evolution. 

In both post-core-bounce supernovae and an early 
universe with net lepton numbers, local net neutrino 
number densities can exceed electron and baryon num- 
ber densities. These large neutrino number densities or 
fluxes, sometimes referred to as a "neutrino background," 
require that the usual Mikheyev-Smirnov-Wolfenstein 
(MSW) formalism |lj for calculating the evolution of 
neutrino flavors be modified to include the effects of 
neutrino-neutrino forward scattering. Though the result- 
ing problem of neutrino flavor evolution can be compli- 
cated, we can identify a key parameter governing the rel- 
evant physics: the ratio of the product of neutrino energy 
and the off-diagonal potential from neutrino-neutrino for- 
ward scattering relative to the difference of the squares of 
the vacuum neutrino mass eigenvalues. If this ratio be- 
comes very large as a result of adiabatic neutrino flavor 
evolution, we can even find an interesting solution, which 
has neutrinos and antineutrinos simultaneously maxi- 
mally mixed in medium. This is different from the usual 
MSW case where, at a given location or time, neutrino 
flavor mixing in medium is maximal for a specific range of 



neutrino energies (this range is narrow for small vacuum 
mixing angles), while antineutrino mixing is suppressed; 
or vice versa. As we discuss below, the off-diagonal po- 
tential from neutrino-neutrino forward scattering plays a 
unique role: it can alter neutrino flavor evolution into a 
form that is utterly unlike the MSW case. 

In the presence of a significant neutrino background, 
neutrino flavor histories can be followed by solving a 
mean-field Schrodinger-like equation in the modified- 
MSW format. This is a long-standing and vexing prob- 
lem. In fact, it has defied general and complete solution 
for the supernova environment, even with sophisticated 
numerical treatments. The existence and importance 
of the flavor-diagonal potential from neutrino-neutrino 
forward scattering and how it might modify MSW-like 
neutrino flavor evolution in supernovae were pointed out 
early on [2] (see Ref. Q for a subsequent formal treat- 
ment). However, the existence of the corresponding off- 
diagonal potential was established only later p|. This 
latter discovery may prove to be a watershed event in 
supernova neutrino physics. 

There have been several attempts to elucidate 
how flavor-diagonal and/or off-diagonal potentials from 
neutrino-neutrino forward scattering can affect active- 
active neutrino flavor transformation in the post-core- 
bounce supernova regime, especially as regards shock re- 
heating [E1S0 and r-process nucleosynthesis 
The rationale for these studies was that the energy spec- 
tra and/or the fluxes of the various neutrino flavors could 
differ on emergence from the neutron star surface or neu- 
trino sphere, and therefore, flavor inter-conversion above 
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this surface could alter these spectra and/or fluxes to 
change supernova dynamics and nucleosynthesis or the 
neutrino signal in a detector. If the energy distribution 
functions and the associated net energy luminosities are 
the same for all flavors of neutrinos and antincutrinos, 
then, obviously, flavor transformation will have no effect. 
However, on account of core de-leptonization and con- 
comitant changes in composition, size, and equation of 
state, this is unlikely to be the case over the entire post- 
core-bounce period of ~ 20 s during which neutrino fluxes 
are appreciable. 

If at any epoch in the supernova environment there de- 
velops a hierarchy of average neutrino energies or lumi- 
nosities among the different neutrino flavors, then flavor 
conversion could alter the rates of electron neutrino and 
antineutrino capture on free nucleons: 

v e + n — » p + e~, (1) 
v e +p -> n + e + . (2) 

These are the processes principally responsible for de- 
positing energy in the material behind the shock after 
core bounce. Therefore, altering their rates by, for ex- 
ample, swapping flavor labels between possibly less en- 
ergetic electron neutrinos and more energetic mu and/or 
tau neutrinos could significantly affect the prospects for 
a supernova explosion Also, the competition be- 

tween these processes and their reverse reactions sets the 
neutron-to-proton ratio in neutrino- heated material y|. 
In turn, this ratio is sometimes a crucial parameter for 
r-process 0, @ and other heavy-element nucleosynthesis 
[ll| associated with slow neutrino-heated outflows. 

Most of the studies cited above posited the existence 
of neutrino mass-squared differences > 0.2 eV 2 . This was 
required for normal MSW resonances to occur in the 
high-density regions most relevant for supernova shock 
re-heating and r-process nucleosynthesis. These regions 
lie above but relatively close to the neutron star, gener- 
ally within a few hundred kilometers. Without the hy- 
pothesized high neutrino mass-squared differences, con- 
ventional MSW evolution in these regions would not re- 
sult in any significant neutrino flavor conversion. 

Although we do not know the absolute vacuum neu- 
trino mass eigenvalues, mi, m^, and ma, the two inde- 
pendent differences of their squares are now measured to 
be Sm 2 w 7 x I0~ 5 and 3 x 1 0~ 3 eV 2 by observations 
of solar and atmospheric neutrinos, respectively. The 
lower 8m 2 has also been measured directly by the Kam- 
LAND reactor experiment. (See Ref. 0| f° r a review of 
neutrino properties.) As these Sm 2 values are certainly 
small compared with the scale previously believed to be 
most relevant for supernovae, one may tend to conclude 
that they have no consequence for supernova dynamics 
and nucleosynthesis. However, as we will discuss below, 
neutrino background effects could alter this conclusion 
dramatically. 

Similar to the supernova case, when there are net 
lepton numbers residing in the neutrino sector in the 



early universe, neutrino flavor conversion can be impor- 
tant in, for example, setting the neutron-to-proton ratio 
and, hence, the 4 He abundance yield in primordial nu- 
cleosynthesis. It was recognized in Ref. that the 
flavor-diagonal potential from neutrino-neutrino forward 
scattering affects neutrino propagation in the coherent 
limit of the problem. However, for reasons that will be- 
come clear, a complete treatment of active-active neu- 
trino flavor conversion in the early universe requires a 
coupled calculation including both flavor-diagonal and 
off-diagonal potentials from neutrino-neutrino forward 
scattering as well as general inelastic neutrino scattering 

@ll5l ll6l| . In fact, the seminal numerical work in Refs. 
Hq shows that the measured neutrino mass-squared 
differences and mixing angles result in an "evening up" 
of the initially disparate lepton numbers residing in neu- 
trinos of different flavors. The flavor oscillations in these 
calculations exhibit "synchronization" in time and space 
and can correspond to near-maximal flavor mixing for 
neutrinos and antineutrinos. We will argue below that 
this numerical result is closely related to our solution 
for the special case where the off-diagonal potential from 
neutrino-neutrino forward scattering dominates. 

In Sec. II we will outline how neutrino flavor trans- 
formation proceeds in the coherent limit when neutrino 
backgrounds are non-negligible. In Sec. Ill we discuss the 
particular limit of domination by the off-diagonal poten- 
tial from neutrino-neutrino forward scattering and the 
corresponding solution. We also speculate under what 
conditions and to what extent this solution could be at- 
tained. Similar issues for a lepton-degenerate early uni- 
verse are discussed in Sec. IV. Conclusions are given in 
Sec. V. 



II. COHERENT FLAVOR EVOLUTION WITH 
NEUTRINO BACKGROUNDS 

Here we give a brief synopsis of coherent neutrino fla- 
vor amplitude development in the supernova and early 
universe environments. In the supernova core and the 
dense environment immediately above it, and in the 
early universe prior to weak decoupling, non-forward neu- 
trino scattering can result in neutrino flavor conversion 
through de-coherence. We will ignore this in what follows 
and instead concentrate on the purely coherent evolution 
of the neutrino fields. It should always be kept in mind, 
however, that our considerations may need to be modi- 
fied at high density or in high neutrino flux regimes. 

Even in the purely coherent limit, following the effects 
of neutrino-neutrino forward scattering in the most gen- 
eral case is daunting in scope. Part of the difficulty in 
following neutrino flavor evolution in the supernova envi- 
ronment is geometric: flavor evolution histories on differ- 
ent neutrino trajectories are coupled. This is because two 
neutrino states will experience quantum entanglement to 
the future of a forward scattering event occurring at the 
intersection of their world lines. In light of this entan- 
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glement there has been considerable speculation about 
whether neutrino flavor evolution can be modeled ade- 
quately b y a mean-field treatment with Schrodinger equa- 
tions 17, 181. The mean field in this case is the potential 
seen by a neutrino by virtue of forward scattering on par- 
ticles in the environment that carry weak charge. Here 
we will follow the conclusions of Ref. 01 an d take the 
mean-field treatment as sufficient in a statistical sense. 
This seems reasonable for the early universe and super- 
nova environments because a statistically large number 
of neutrino scattering events and entanglements occur in 
these places. 

Further complicating the supernova problem is the 
non-isotropic nature of the neutrino fields above the 
neutrino sphere. Neutrinos traveling along trajectories 
nearly tangential to the neutrino sphere may have quite 
different flavor amplitude histories from those moving ra- 
dially or near radially. For now we will ignore this feature 
of flavor development and instead approximate all neu- 
trinos as evolving the way radially propagating neutrinos 
evolve. This approximation is not a good one (it has been 
made in all previous numerical work j(| Il0| ) , but it will 
suffice in our analytic arguments here. 



A. Overview 

In vacuum, the flavor (weak interaction) eigenstates of 
neutrinos are related to the mass (energy) eigenstates by 
a unitary transformation U m : 




= U„ 




(3) 



where the mass eigenstates \i>i), \v2), and correspond 
to the vacuum mass eigenvalues, mi, mi, and 7713, respec- 
tively. The unitary transformation U m can be written in 
terms of a sequence of rotations, 

U m = U 23 U 13 U 12 . (4) 

A convenient representation for these rotations is 





u- 



1:; 



Uv 



1 

cos 6*23 sin 6*23 
-sin 6*23 cos 6*23 
cos 9i 3 e lS sin #13 

10 
„-«o;~^„ cos 6»i3 

sin 812 

COS 012 

1 



(5) 




In the above representation, the mixing angles, 612 and 
623, have been measured by observations of solar and at- 
mospheric neutrinos and related experiments. In partic- 
ular, the best fit for the vacuum mixing of the mu and tau 
neutrinos is very near maximal, which gives 823 ~ tt/4. 



However, the mixing angle, #13, and the CP-violating 
phase, 5, have not been measured yet. 

Here we consider a neutrino mixing scenario where 
m 3 > m 2 > mi with 8m\ 2 = m\ — m\ w 7 x 10~ 5 eV 2 
and Sml 3 = m\ - m\ 3 x lCT 3 eV 2 , 6» 23 = 7r/4, and 
5 = 0. With the definitions 



K) = ^ 



V2 ' 

Vu) + Wr) 



V2 

it is straightforward to show that 
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(6) 
(7) 

(8) 



where for example, C12 = cos #12 and s\ 2 = sin #12. The 
states \v*) and \u*) are still useful in medium. This is 
because in the supernova medium the mu and tau neutri- 
nos have very nearly the same interactions, so that mat- 
ter effects on mixing and effective mass for these species 
are nearly identical (likewise for the mu and tau antincu- 
trinos). This will also be true for the early universe if 
the net muon and tau lepton numbers are identical. For 
the sake of our arguments here, we will take the sym- 
metry between mu and tau neutrinos and that between 
their antiparticles to be rigorously true so that \u*) and 
\v*) are effective flavor eigenstates in medium. Since 
6m 2 3 3> <5to 2 2 , the regions of neutrino flavor mixing gov- 
erned by these parameters should be well separated. As 
neutrinos propagate outward from the neutrino sphere 
at high density in supernovae, Sm 2 3 «3x 10~ 3 eV 2 be- 
comes relevant first. We will focus on neutrino flavor 
mixing with this parameter, for which v* is effectively 
decoupled [see Eq. JSj] and we only need consider mixing 
of v e and v* [iji |2Cj. Thus, the general problem of iv 
mixing in medium is reduced to one of 2v mixing in our 
scenario. 

With the above simplification we can hereafter follow 
the notation of Ref. gg. In particular, we now simply 
refer to \v*) as \v T ) and write the effective 1v unitary 
transformation in vacuum as 



\v e ) — cosv\v\) + sintf|i/ 2 ; 
\v T ) = — sa\0\vi) + cos #^2) 



(9) 
(10) 



where and |i/ 2 ) refer generically to the light and 
heavy mass eigenstates, respectively, and 6 is the ef- 
fective 2v vacuum mixing angle. The relevant vacuum 
mass-squared difference is Sm 2 w3x 10~ 3 eV 2 . The cor- 
responding effective vacuum mixing angle is ~ #13 , with 
the current reactor experiment limit being sin 2 2#i3 < 0.1 
(see e.g., Ref. H3). 

Consider a neutrino of initial flavor a — e or r. As it 
propagates outward from the neutrino sphere in super- 
novae, the evolution of its state, the ket (t)}, can be 
described as 



|tf Va (t)) =0tea(i)K)+OW(i)K) 



(11) 
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where a ea (t) or a Ta (t) is the amplitude for the neutrino to 
be a v e or v T , respectively, at time t. (Note that t could be 
any Affine parameter such as radius along the neutrino's 
world line.) Alternatively, the evolution of \^u a {t)) can 
be described as 

|*v„(t)) = a la {t)\ Vl (t)) + a 2a (t)\v 2 {t)), (12) 

where and \v 2 {t)) are the instantaneous mass (en- 

ergy) eigenstates in medium, and ai a (t) and a 2a (t) are 
the corresponding amplitudes. The flavor eigenstates are 
related to |^i(t)) and \v 2 (t)) as 



cos0 M (t)M*)> +sm0 M (t)\v2(t)) (13) 
-sm6 M (t)\Mt)) +cosO M (t)\v 2 (t)), (14) 



where 0M(t) is the effective 2v mixing angle in medium 
at time t. In matrix form, the ket \^i/ a {t)} can be repre- 
sented by 



«(*) 
«(*) 



(15) 



in the flavor basis and by 



a la (t) 

d2a(t) 



(16) 



in the energy basis. In analogous fashion we will employ 
a 2D scheme to follow separately the flavor evolution of 
the antineutrino sector. 



B. Characterizing Neutrino Densities 



The single neutrino density operator at time t pro- 
jected into the energy basis is 



l^(*)K^«(i)l 



K«WI>iW>MOI + \a2 a {t)\>2{t)){v 2 {t)\ 

oia (t)a* 2a (t) I v x (*)) {v 2 (t) | + a* Xa (t)a 2a (t) \ v 2 {€)) {v x (t) \ . 

I 



(17) 



The second line of Eq. (|17|l contains cross terms which in 
general have complex coefficients. However, these cross 
terms vanish in the limit where neutrino flavor evolution 
is adiabatic. This is because a neutrino evolving adiabat- 
ically is always in a single energy state. For example, in 
this limit we might have |ai Q (i)| = 1, which would imply 
that |a2a(i)| = due to the normalization condition 



<* a (t)|* Q (t)) = \a la (t)\ 2 + \a 2a (t)\ 2 = 1. 



(18) 



The density operator for the neutrinos or antineutrinos 
with momentum centered around p in a pencil of neu- 
trino or antineutrino momenta and directions d 3 p can be 
defined as in Ref. |(|: 

P P (t)d 3 P = X^J^WK^WI, (is) 

a 

p p (t)d 3 p = J2 dn ^P a (t))(^P a (t)\- (20) 

a 

Note that the traces of these operators over neutrino fla- 
vor do not give unity but rather the total number density 
of neutrinos or antineutrinos of all kinds in the pencil. 

We assume that neutrinos and antineutrinos of all fla- 
vors are emitted from the same sharp neutrino sphere of 
radius R v in supernovae. (This is not a particularly good 
approximation for neutrinos very near the neutron star 
surface, but it will suffice for our arguments.) At a radius 



r > R v , the neutrino sphere subtends a solid angle of 

AO„(r) = 27r (i-y/i-Rl/r*). (21) 

Within this solid angle, the number density of v a in a 
pencil of directions and momenta is 



dn v „ = 



1 



where L Ua is the energy luminosity of u a , d£l u is the 
pencil of directions, E v is the neutrino energy, f Va (E v ) is 
the normalized energy distribution function for v a , and 
(E Ua ) is the corresponding average v a energy. Here and 
in the rest of this paper we assume that neutrinos have 
relativistic kinematics and employ natural units where 
ti = c = 1. The function f Va {E v ) can be fitted to the 
results from supernova neutrino transport calculations 
and is commonly taken to be of the form 



fu a {E v ) 



El 



1 



(23) 



where T Va and r\ Va are two fitting parameters and F 2 {rj Va ) 
is the Fermi integral of order 2 and argument r\ Va . The 
Fermi integral of order k and argument T) is defined as 



<: k dx 



e x- V + I 



(24) 
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In terms of these integrals, the average v a energy is 



poo 

(E v<x ) = / E v f„ a (E v )dE v =T Ua 
Jo 



(25) 



We assume that all neutrino species have thermal, 
Fermi-Dirac energy distribution functions in the early 
universe, so the number density of v a in a pencil of di- 
rections and momenta is 



dn v 



1 [dn„\ E 2 dE u 



2tt 2 V 4tt J e E "/ T "c l -v« a +l' 



(26) 



Note that although the above equation uses the same 
symbols T„ Q and r\ Va as Eq. (|2*S|l . the physical meanings 
of these symbols are very different. In Eq. (|23[) . T„ Q and 
rju a are simple parameters used to fit the energy distribu- 
tion functions obtained from supernova neutrino trans- 
port calculations, while in Eq. (|26() . T Va is the tempera- 
ture and r\ Va is the degeneracy parameter of the v a gas 
in the early universe. However, as the assumed neutrino 
energy distribution functions for the supernova and early 
universe environments have the same functional form, we 
use the same symbols in both cases for convenience. For 
the homogeneous and isotropic neutrino gas in the early 
universe, the local proper number density of v a is 



2^3 
2ir z 



(27) 



Equation l|26|) can be rewritten as 



dn,j 



47T 



/„„ (E v )dE„. 



(28) 



In terms of the scaled neutrino energy e = E v jT Va ^ 



f Va {E v )dE v = f Va (e)de 



F 2 {r]u a ) +l 



(29) 



In analogy to the baryon-to-photon ratio r\ = 
(rib — Tii) /n 7 «6x 10 -10 , we define a z^-to-photon ratio 



n,, - til 



ii 



7T 2 (T Va 

12C(3)VT 7 



(30) 

where T 7 is the photon temperature, ( (3) = 1.20206, and 
we have used T Pq = T Va and r\v a = —r] Va to obtain the 
last identity. For T„ a — T 7 and small £ Va , t] Vot rs 1.46£v a . 
Current limits on all lepton numbers are £ Va < 0.1 (cf. 
Ref. 0). 



C. Neutrino Propagation in Medium 

For a neutrino originating as a v a at t — 0, its subse- 
quent flavor evolution along a radially-directed trajectory 
with Affine parameter t is described by 



|*0 = (#vac + H ev + ff„ v ) (31) 



where we have decomposed the overall evolution Hamilto- 
nian into contributions from vacuum neutrino masses and 
from mean-field ensemble averages for neutrino-electron 
and neutrino-neutrino forward scattering. These contri- 
butions are discussed individually below. 

For a neutrino with energy E v and vacuum mass m <C 
E v , we have E v = ^/p 2 + m 2 w p + m 2 /(2p), where 
p is the magnitude of the neutrino momentum p. In 
this limit, the vacuum-mass contribution to the flavor 
evolution Hamiltonian is 



pi 



1 

Yp 



'm 2 \^)(^\+ m 2 W 2 )(^\), (32) 



where / is the identity operator. 

Electron neutrinos and antineutrinos can forward scat- 
ter on electrons and positrons through exchange of W ± . 
In contrast, there is no such charged-current forward 
scattering for v^, u T , and their antiparticles due to the ab- 
sence of /i^ and in the environments of interest here. 
Consequently, the effective contribution from charged- 
current neutrino-electron forward scattering to the flavor 
evolution Hamiltonian is 



H w (t) = A(t)\v e )(v e \, (33) 



where 



A(t) = V2G F {n e - -n e +) = V2G F n b Y e . (34) 

In the above equation, ri e - , n e + , and are the proper 
number densities of electrons, positrons, and baryons, re- 
spectively, at the position corresponding to time t, and 
Y e = [n e - — n e +) /rib is the net electron number per 
baryon, or electron fraction. 

For a specific neutrino with momentum p, the effec- 
tive neutral-current neutrino-neutrino forward scattering 
contribution |2( to the flavor evolution Hamiltonian is 

H vv (t) = V2G F J {I- cos0 pq ) [p q (i) - q (t)] d 3 q, 

where q is the momentum of the background neutrinos 
and cos # pq = p-q/pq. The term (1 — cos # pq ) stems from 
the structure of the weak current 2] . This can be seen 
from the limit where completely relativistic neutrinos are 
traveling in the same direction along the same spacetime 
path. In this limit 1 — cos # pq = and neutrinos never 
forward scatter on one another. Obviously, for the ho- 
mogeneous and isotropic neutrino distribution functions 
characteristic of the early universe, cos# pq averages to 
zero and the ensemble average of (1 — cos# pq ) is unity. 
In the supernova environment the term (1 — cos # pq ) will 
be largest close to the neutron star, where the neutrino 
trajectories can intersect at high angles. At sufficiently 
large radii above the neutron star, the neutrino-neutrino 
forward-scattering contribution to the flavor evolution 
Hamiltonian will scale as r~ 4 . As the neutrino-electron 
forward-scattering contribution will scale roughly as r~ 3 , 
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it may be dominated by the neutrino-neutrino forward- In matrix form, the neutrino flavor evolution equation 
scattering contribution at small to moderate distances in the flavor basis is 
from the neutron star. 



0*/ 

'"dt 



p- 



4p 



A 

2 +<*> 



- 1 
/+ 2 



A + B-Acos26> Asm20 + B er 
A sin 29 + B Te A cos 29 - A - B 



(36) 



where we have separated the Hamiltonian into a traceless term and a term proportional to the identity matrix. The 
latter term gives only an overall phase to the neutrino states, and is therefore unimportant in neutrino flavor conversion. 
In the above equation, A = 5m 2 /2E l/ , and a„, B, and B eT (B Te — B\ r ) are the potentials from neutrino-neutrino 
forward scattering. Specifically, 



V2 



GfJ (1 - cos0 pq ) ([p q (t) - ^ q (t)] ee + [p q (i) - ^ q (i)] rr )d 3 q, 



B = V2G F 
B er = 2\/2Gf 



(i 



cos I 



(1 



q) [PqW-P q W] er d 3 q, 



where the matrix elements of the density operators are defined as 

[p q (t)-p q (t)] TT rf 3 q = H,3 q (t)d 3 qK) 

[pq(«)-p q (*)] eT ^ 3 q = (^e|/3 q (t)d 3 q|^) 

I 



(P e |^ q (t)d 3 q|P e ), 
(P T |/5 q (£)d 3 q|P r ), 

(^e|p q Wd 3 q|^r). 



(37) 
(38) 
(39) 



(40) 
(41) 
(42) 



The physical interpretation of these matrix elements is 
straightforward, even if the notation is cumbersome. 
For example, [pq(t) — Pq(t)] ee ^ 3 q gives the expectation 
value for the net v e number density in the pencil of mo- 
menta and directions cZ 3 q centered on q. Note that the 
off-diagonal matrix element vanishes and makes no con- 
tribution to B eT if neutrinos remain in their initial flavor 
states. This is evident if we expand out the first term in 
Eq. O: 



(43) 



In the above equation, one or the other amplitude in 
the sum on the right-hand side will be zero unless some 
neutrino flavor transformation has occurred at the time 
t when this matrix element is evaluated. 

For real B eT = B Te , it is convenient to define the ef- 
fective mixing angle 9m in medium by 



cos 26 M (t) 
sm29 M (t) 



(Acos26> - A 
(A sin 26» + B e 



- B)/A eS , 
)/A cff , 



(44) 
(45) 



where 



A cff = \J (A cos 26» - A - B) 2 + (A sin 29 + B eT f . (46) 
With the term proportional to the identity matrix 



dropped, Eq. I|36|l can be transformed to the instanta 
neous energy basis to give 

, 9^m = [ -A off /2 -i6 M (t) 
dt [ i9 M (i) A eff /2 



(47) 



where 9m (i) = (19m /dt. In the limit where \9 m (t) | <C 
Aeff/2, Eq. (|4TI) becomes two decoupled equations and 
flavor amplitude evolution is adiabatic. The flavor evo- 
lution equations and the corresponding effective mixing 
angle 9m in medium for antineutrinos can be obtained 
from those for neutrinos by replacing A, B, and B eT in 
the latter with —A, —B, and —B eT , respectively. 

The condition \9 m (t) | <C A er r/2 for adiabatic neutrino 
flavor evolution is most stringent when A e g reaches the 
minimum value | A sin 29 + B eT | at an MSW resonance 
corresponding to 

A cos 26 1 = ,4 + 5. (48) 

At resonance, the effective in-medium mixing angle is 
(ires) = tt/4 and mixing is maximal with sin 2 26 a 



'M 



1M 



1. We can define an adiabaticity parameter 
A cff (i rcs ) (Asin26» + B eT ) : 



where 



2\9m (ti 



H 



Acos26> 



-n, 



V 




A - 


- B 


V 


res 


Ah 


- B 



(49) 



(50) 
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is the scale height for the total potential V = A + B at 
resonance with V — dV/dt. We can gain more insight 
into the adiabaticity parameter by further defining a res- 
onance region corresponding to 1/2 < sin 2 29m < 1. In 
this region the change in V around the resonance value 
A cos 20 is 5V = | A sin 20 + B eT \, so the width of this 
region is 



SV |Asin20 + 5 e 



A cos 29 



-H. 



As the oscillation length at resonance is 
2tt 2tt 



A e ff (t rc 



|Asin20 + 5 e 



we have 



7 = 2tt 



(St); 



(51) 



(52) 



(53) 



In summary, large B eT increases 7 in two ways: (1) by 



increasing the resonance width (5t) Tes and; (2) by de- 
creasing the oscillation length L rcs . 

Clearly, neutrino flavor evolution will be adiabatic for 
7 ^> 1. For the small effective vacuum mixing angle 
9 <C 1 of interest here, neutrino flavor conversion will 
be complete in this limit. For arbitrary 7, the prob- 
ability of neutrino flavor conversion after propagation 
through resonance is well approximated by 1 — Plz, where 
Plz = cxp (— 7T7/2) is the Landau-Zener probability for a 
neutrino to jump from one energy eigenstate to the other 
in traversing the resonance region. 



D. Neutrino Potentials in the Adiabatic Limit 

Due to the cross terms in the single neutrino density 
operator in Eq. (|17fl . B eT is generally complex. If these 
cross terms are unimportant, then both B and B eT are 
real and their expressions in Eqs. (|38|) and (|39|l can be 
simplified as 



B 



B e 



= -V2G F J2 / (l-cos0 pq ) 
= V2G F J2 I (l-cos0 pq ) 



cos29m (l — 2|ai Q | 2 ^ dn Va — cos 29m fl — 2|ai Q | 



dn,-, 



sin 29m ( 1 ~ 2|ai a | 2 ) dn v — sin 29m ( 1 — 2|ai Q | 2 ) dn^ 



(54) 
(55) 



As mentioned in Sec. Ilffjl the cross terms in Eq. H17|) vanish if neutrino states evolve adiabatically. In this limit, 
the above expressions of B and B eT can be simplified further. As the electron and neutrino number densities at the 
neutrino sphere are far above those satisfying the resonance condition, v e and v T are born essentially as the energy 
eigenstates \v 2 ) and \v\), respectively. For adiabatic evolution, \ai e \ 2 w and |ai T | 2 w 1 for all subsequent time t. 
For antineutrinos, adiabatic evolution gives \a\ e I pa 1 and I a iT I ~ 0. Thus, we have 



B pb V2G F J (l-cos0 pq ) [{dn VT 
B eT « V2G F J (I — cos P q) \{dn Vl , 



dn Vll ) cos 29m + [dnp T — dn^^ ) cos 20a/] , 
dn VT ) sin 20m + (dnp e — dri[, T ) sin 29m] ■ 
I 



(56) 
(57) 



III. FLAVOR MIXING WITH LARGE 
OFF-DIAGONAL POTENTIAL 

We now focus on adiabatic neutrino flavor evolution, 
for which the potentials from neutrino-neutrino forward 
scattering are given by Eqs. (|56|l and Q57J1. Our main 
concern is the effects of these potentials on neutrino fla- 
vor evolution in supernovae. In Sec. Ill Cl we have outlined 
this evolution for a specific neutrino as it propagates out- 



J 



ward from the neutrino sphere. At a given radius r with 
potentials A and B, the resonance condition in Eq. I|48() 
will be met for a particular neutrino energy £V CS , i.e., 



cos 20 — A + B. 



2E n 



(58) 

As B and B eT at a given radius involve integration of 
cos 20 Mi cos 29 Mi sin20 M , and sin20 M over the neutrino 
energy distribution functions, it is convenient to write 



cos 29 M — 



1 — E v J E T , 



y/(l - E„/E rcs y + [tan 20 + (E v /E rcs )(2E YCS B eT )/(Sm 2 cos20)p 



(59) 
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tan 20 + (E v /E, cs )(2E, cs B er )/(5m 2 cos 20) 



sin 29m = 
cos 29m — 
sin 29m = 





- E v /E Tes f - 


- [tan 20 + 


- {E u /E rcs ){2E rcs B eT )/(8m 2 


cos20)] 2 






1 - 


f E„/E ICS 






f- E v /E ies ) 2 - 


- [tan 20 - 


(E v /E ICS )(2E rcs B eT )/{Sm 2 


cos20)] 2 




tan 20 - 


- (E v /E Tei 


i )(2E ICS B eT )/(5m 2 cos 20) 






f- E v /E ves f - 


- [tan 20 - 


(£,/£ rcs )(2£; rcSJ B er )/(c5m 2 


cos20)] 2 



r 



Let us now consider the limits of the above expressions 
for the in-medium mixing angles when 9^1 and B eT is 
positive and so large that the second terms in the square 
root of these expressions dominate the first terms. In this 
limit we have 



<C 1. Furthermore, B eT also affects the in- medium 
mixing for antineutrinos, which is strongly suppressed 
(sin 2 20m <C 1) in the absence of neutrino-neutrino 
forward scattering. For B eT = (5m 2 /2E YCS ) cos 20 and 
<C 1, antineutrinos with E v > E les have substantial 
in-medium mixing with 1/5 < sin 2 20m < 1/2. 



cos 20m - 


+ o, 


(63) 




sin 20m - 


- 1, 


(64) 




cos 29 m ~ 


- o, 


(65) 


A. Towards a Self-Consistent Solution with a 


sin 20 A / - 


- -1, 


(66) 



for which both neutrinos and antineutrinos have maximal 
in-medium mixing with 



77 

4' 
3tt 

T' 



(67) 
(68) 



A large negative B eT clearly changes the signs of the lim- 
its in Eqs. (|63fl - l|66|) . and the in-medium mixing angles 
in this case are 0m — > 37r/4 and 0m — > 7r/4. 

In general, we see that large in-medium mixing will 
occur simultaneously for neutrinos and antineutrinos over 
a broad range of energies if adiabatic flavor evolution 
results in 



|£ eT | » 



2E r , 



■ cos 20 



(69) 



at some radius above the neutrino sphere. Using Eq. (1581) . 
we can rewrite the above equation as \B eT \ ^> A + B. As 
cos20m — * and cos 20m — > when this is achieved, 
Eq. gives B — > 0. Therefore, Eq. ||BSJ| reduces to 



I Berl > A 



(70) 



Note that even a | fieri on ly as large as 
(dm 2 /2E les ) cos 20 already has important effects on 
the in-medium mixing of neutrinos and antineutrinos. 
While a neutrino with resonance energy E v = E Tes has 
maximal in- medium mixing independent of B eT , the 
energy range over which neutrinos have large in-medium 
mixing with 1/2 < sin 2 20m < 1 is strongly affected 
by B eT . For B eT = 0, this energy range corresponds 
to E Tes (l - tan 20) < E u < E IeB (l + tan 20), which 
is very narrow for J < 1. In contrast, for example, 
with B eT = (Sm 2 /2E res ) cos 20, all neutrinos with 
E v > E ICS /2 have 1/2 < sin 2 20 M < 1 even for 



Here we outline a possible self consistent B er -dominant 
solution (BDS) which meets two criteria: (1) |fi eT | ^S> A; 
and (2) adiabaticity, 7 > 1. An immediate question is: 
can adiabatic neutrino flavor evolution ever produce a 
large B eT as in Eq. (|7U|) ? Obviously, the answer is yes 
if one can demonstrate that this result is obtained under 
some conditions. Such demonstration requires following 
the flavor evolution of neutrinos with a wide range of 
energies covered by their energy distributions. As men- 
tioned in Sec.[n] this process may sound straightforward 
but turns out to be computationally difficult. On the 
other hand, if adiabatic flavor evolution can indeed give 
rise to a large B eT , then the conditions required for this to 
occur must depend on the following: the neutrino mixing 
parameters 5m 2 and sin 2 20, the profile of electron num- 
ber density that gives the potential A, and the neutrino 
luminosities and energy distribution functions that are 
related to the potentials B and B eT . Our goal here is to 
examine these dependences. In so doing, we will not be 
able to answer the question posed at the beginning of this 
paragraph, but we will be able to provide a range of con- 
ditions that can guide future numerical calculations in 
search of a complete solution for neutrino flavor mixing. 

We start with the basic input for our discussion. As 
explained in Sec. Ill Al the mixing parameters of interest 
here are Sm 2 « 3 x 10" 3 eV 2 and sin 2 20 < 0.1. To 
characterize the potential A, we need the electron number 
density n e — Y e nt,. We note that the envelope above the 
post-core-bounce neutron star can be approximated as 
a quasi-static configuration with a constant entropy per 
baryon S in the gravitational field of the neutron star. 
In this case, the enthalpy per baryon, TS, is roughly 
the gravitational binding energy of a baryon, so that the 
temperature T scales with radius as 



rp _ M NS m p s -i r -i 



m ■ 



(71) 



pi 
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where mp\ ~ 1.221 x 10 22 MeV is the Planck mass, m p is 
the proton mass, and Mns is the neutron star mass. At 
late times relevant for r-process nucleosynthesis, the envi- 
ronment above the neutron star is radiation-dominated, 
so 



n 2tt 2 T 3 
S « — g s — 
45 rib 



(72) 



in units of Boltzmann constant kg per baryon. In the 
above equation, g s is the statistical weight in relativistic 
particles: g s ~ 11/2 when e^-pairs are abundant and 
g s w 2 otherwise. Combining Eqs. (|71|l and l|72|) . we 
obtain the run of baryon number density for the r-process 
epoch as 



27r 2 / M N sm p , 
45 V m F 



'pi 



(73) 



The potential A is given by 
2 v / 2tt 2 



rib 



45 



g s Y e G F { 



'pi 



« (5.2xlO- 13 MeV) 5s y e ( I ^) S^ 3 , (74) 

where Sioo is 5 in units of lOOfcs per baryon and rg is r 
in units of 10 6 cm. The r-process epoch corresponds to a 
time post-core-bounce t p b > 3 s. This is a relatively long 
time after core bounce, at least compared with the time 
scale of the shock re-heating epoch and the time scale 
for evolution of neutrino emission characteristics such as 
luminosities and average energies. The potential A in 
Eq. (|?4"|) can also be used to describe crudely the shocked 
regions of the envelope above the core in the shock re- 
heating epoch, £ p b < 1 s, if we take g s ~ 1 and employ a 
low entropy 

To evaluate B and B eT , we need the differential number 
density of each neutrino species at radius r > R v above 
the neutrino sphere. The differential v a number density 
in the absence of flavor evolution is given by Eq. (|22H . 
which depends on the luminosity L Va and the average 
energy (E V J). For some illustrative numerical estimates 
we will assume that all neutrino species have the same 
luminosity, 



(75) 



and take the average neutrino energies to be 

(E Ue ) = 10 MeV, {E Be ) = 15 MeV, 

{E Ut ) = (E„ r ) = 27 MeV. (76) 

Assuming that adiabatic flavor evolution up to some 
radius r > R v results in a BDS described by the cri- 
terion in Eq. JTUJ), we now examine the implications of 
this criterion for supernova conditions. For dcfmitcncss, 
we discuss the case of a large positive i? B r DS , which gives 
sin 29m — * 1 and sin26>M — ► — 1. In this case, Eq. i|5*7|) 
gives 



B L B T DS « V2G F J (1 - cos6» pq ) (dn Vr - dn^J, 



(77) 



where we have assumed that the luminosities and en- 
ergy distribution functions for v T and v r are very nearly 
the same in the supernova environment. This is a good 
approximation because these species experience nearly 
identical interactions both in the dense environment of 
the core and in the more tenuous outer regions. For a 
radially-propagating test neutrino, the intersecting an- 
gles of the background neutrinos, # pq , are coincident with 
the polar angle in the integration over d£l v for the test 
neutrino. Assuming that neutrinos of all flavors originate 
on the same neutrino sphere and using dn t/c and dn Pi! of 
the form in Eq. 5f T DS in Eq. O can be evaluated 



B 



BDS 



V2G F 
A-kRI 



1 - y/1 - Rjp 



(2.1 x 1(T 10 MeV) R-l (\ - ^l-i? 2 6 /r 2 ) 



L 



L 



p, o2 



(E Ue )/(WMeV) (E Sa )/(10MeV) 



(78) 



where R v6 = i?„/(10 6 cm), L Ue52 = ^/(lO 52 ergs s" 1 ), 
and L Pe 5 2 = L Se /(10 52 ergs s" 1 ). For r > R v ~ 10 6 cm 
and the assumptions in Eqs. I|75|) and J7BJ, we obtain 



B 



BDS « (1.8 x 10- 11 MeV) ^ 52 ^ 6 . (79) 



Using Eqs. I|74|) and (|78|l . we can rewrite the criterion 
B BDS > A as 



Rl/r 



1 - Vi - WF 2 



45 / to 2 ,) 



(398) 



1.4M 







M NS 



R V S 

9sY e 



RuqS^ 



100 



9sY e 



L 



u c 52 



//, 52 



(-E„ e )/(10MeV) (£ Pe )/(10MeV)_ 



(80) 
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For r ^> R v and the assumptions in Eqs. (|75|) and 1)76(1 . 
we obtain 



r 6 < 33 



1.4M 







Mns 



S W0 L„ 52 R 



3s 



(81) 



We have assumed adiabatic neutrino flavor evolution 
in the above discussion. Though the adiabaticity of the 
general flavor evolution can be ascertained only with a 
sophisticated numerical treatment, the BDS clearly will 
not be self-consistent if flavor evolution is not adiabatic 
at the radius where i?f T DS ^> A is achieved. On the other 
hand, if we can show that at this radius the adiabatic- 
ity parameter for the neutrino with energy E rcs satisfies 
7 BDS ^S> 1, then the BDS is more likely to be obtained. 
With H \A/A\ w r/3, this criterion [see Eq. can 
be rewritten as 



„BDS 



B ernj 

15 



'PI 



M NS m p 
( L ^ 



1 - y/1 - Rl/r 2 
Ru/r 



g s Y e J \{E Ve ) (E Pe ) 



> 1, 



(82) 



which reduces to 



„BDS 



10' 



1.4M 



NS 



C4 t2 
J 100-^^52 

9sY e 



> 1 (83) 



for r 3> R v and the assumptions in Eqs. 1)75(1 and ((76(1 . 

The criteria in Eqs. 1(801) and ((82(1 can be met in some 
regions with significant scale and duration above the 
neutron star during both the r-process and shock re- 
heating epochs. For example, an r-process environment 
with modest entropy might have Y e w 0.4, R v § w 1, 
g s w 11/2, ^ioo ~ 1-5, and L„52 ~ 0.1. For these pa- 
rameters Eqs. (|5U|) and would give r 6 «8 and 22, 
respectively, so the BDS may be obtained over an ex- 
tended region above the neutrino sphere. For a higher 
entropy, 5ioo = 2.5, but with all the other parameters re- 
maining the same, Eqs. 1(80(1 and 1(82(1 would give r$ <C 60 
and 40, respectively. We can put these limits in perspec- 
tive by noting the temperature at which salient events 
or processes occur above the neutron star. The radius 
corresponding to a temperature Tg (measured in units of 
10 9 K) is very roughly r 6 w 22.5/(T 9 5i o) [sec Eq. lf7TJTj. 
Weak freeze-out, where the neutron-to-proton ratio is set, 
occurs at Tg ~ 10. The neutron capture regime in the re- 
process is typically further out, occuring between Tg ss 3 
and Tg » 1. Therefore, the limits on the radius discussed 
above are so generous that maximal in-medium mixing 
for both neutrinos and antineutrinos associated with the 
BDS could affect important weak interaction processes 
in the envelope, the r-process, and the neutrino signal. 

Taking Y e = 0.4, g s = 11/2, M NS = 1.4 M Q , R v = 
10km, L Ve = L De = L v , (E Ve ) = lOMeV, and (E De ) = 
15 MeV, we use Eq. 1(80(1 to calculate the combinations of 
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FIG. 1: Contours of limiting radius (in units of 10 6 cm) be- 
neath which Bf^ s S> A may be obtained. Except for the 
region corresponding to the larger values of the limiting ra- 
dius, these results are generally more stringent than those 
from 7 BDS 3> 1. The chosen range of parameters is meant to 
be characteristic of the r-process epoch. The horizontal axis 
is entropy in units of lOOfes per baryon, while the vertical axis 
is neutrino luminosity L„ in units of 10 52 ergss _1 . Here we 
take Y e = 0.4, g s = 11/2, M N s = 1.4 M s , and R v = 10km. 
We also assume that L Vll — L Vf = L v , (E v<! ) — 10 MeV, 
(Ei) s ) = 15 MeV, and that all mu and tau neutrinos and an- 
tineutrinos have identical luminosities and energy spectra. 



L v and S for which the criterion B^ r DS ^S> A can be met 
below a fixed radius. The results are shown in Fig. \\\a,s 
contours labelled by the limiting radius. Except for the 
region corresponding to the larger values of the limiting 
radius, these results are generally more stringent than 
those from the criterion 7 BDS ^> 1 . Compared with the r- 
process regime at later times, the shock re-heating epoch 
is characterized by much higher neutrino luminosites. In 
general both L Ue and L Pe are ~ 10 52 ergss _1 . Taking 
Y e = 0.35, g s = 1.5, and R v = 40 km (other parameters 
remaining the same as for Fig. [Q, we present in Fig. [21 
the constraints on L v and S for which Bf^ s 3> A can 
be met below various radii during the shock re-heating 
epoch. Based on these results, both shock re-heating and 
the neutrino signal could be affected by maximal neutrino 
flavor mixing |5|, 12 1| if there were a hierarchy of neutrino 
energies at this epoch. 

Note that average v e and v e energies may be quite 
similar during much of the shock re-heating epoch, but 
the luminosities for v e can be significantly larger than 
those for v e . This is especially true for shock break-out 
through the neutrino sphere, the so-called neutronization 
burst. For a time span of ~ 10 ms we could have L u ^2 ~ 
10, while ip e 52 is an order of magnitude smaller. Since 
neutrino flavor mixing in the coherent limit is a phase 
effect, the 10 ms duration of this high-luminosity burst 
may be long enough to establish the BDS. Neglecting 
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FIG. 2: Same as Fig. Q except that now the conditions are 
meant to be representative of the shock re-heating epoch. In 
this case we take Y e = 0.35, g s = 1.5, and i?„ = 40 km. 



L Pe and taking L„ c52 = 10, (E Ve ) = lOMeV, R v6 = 
4, Y e — 0.35, g s = 1.5, and 5ioo = 0.15, we find that 
Bff s > A may be obtained for r 6 < 15 [Eq. JgQJ]. 
This limit becomes rg <C 50 if 6*100 = 0.2. For both 
cases the limit from Eq. l|82[l is much weaker. Therefore, 
the neutrino signal from the neutronization burst and 
the early shock re-heating process could be affected by 
maximal neutrino flavor mixing associated with the BDS. 

B. Is the BDS Ever Attained? 

Achieving the BDS is dependent on a number of con- 
ditions, many of which are unlikely to strictly and gen- 



where we define the effective net number of neutrinos per 
baryon through B = \/2GFribY£ s . Because 5m 2 is small, 
the resonance energy also tends to be small at distances 
where neutrino fluxes are appreciable. 

However, if neutrinos transform their flavors via a 
strict MSW evolution then B (and Y° s ) will drop with 
the radius of the fluid element and eventually will be 
driven negative [|| . To see this consider first an example 
(hierarchical) energy spectrum for v e and v T neutrinos 
as they leave the neutrino sphere. In Fig. @ we show 



erally obtain in environments in nature with high neu- 
trino fluxes. The essence of the BDS is the dominance 
of the flavor off-diagonal potential and, in particular, 
2E v B eT ^> Sm 2 cos 29. Since the measured neutrino 
mass-squared differences are small, it will not take a 
large flavor off-diagonal potential to force the system into 
something like the BDS. 



However, as Eq. (|43[1 shows, a necessary condition for 
B eT to be non-zero at some time/position is that some 
neutrinos must have transformed their flavors there. This 
can be problematic because in both the early universe 
and the post-shock supernova environment a fluid ele- 
ment will evolve from conditions of very high density to- 
ward lower density. For example, the region near the 
neutron star surface is very high density, corresponding 
to high electron degeneracy. This will tend to suppress 
in-medium neutrino mixing. A hydrodynamic flow away 
from the neutron star surface will carry a fluid element 
into regions of lower temperature and density and net 
neutrino fluxes. At large enough radius the neutrino- 
electron potential will scale like A ~ r~ 3 , while the 
flavor-diagonal and off-diagonal neutrino-neutrino poten- 
tials will scale as r~ A . As a result, there may be some 
region where the neutrino-neutrino potentials dominate. 



The neutrino resonance energy experienced in this fluid 
element at radius r will be E Tes = 8m 2 cos29/2(A + B). 
Near the neutron star surface -E res Wlu be extremely 
small. Further out, in an adiabatic and roughly hydro- 
static envelope (notation as in the last section), the res- 
onance energy at radius r will be 



(84) 
(85) 



Fermi-Dirac-type energy spectra for these species, tak- 
ing the neutrino degeneracy parameter for both to be 
r\ Va = 3, and taking average energies (E Ue ) = lOMeV 
and (E Vt ) = 27MeV. The actual supernova neutrino 
energy spectra may differ significantly from these, but 
they serve to illustrate general trends. Note that for our 
chosen spectral parameters, the v e population at lower 
energies is larger than the v T population for comparable 
luminosity in the two neutrino species. 

As our example fluid element moves out to larger r, 
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FIG. 3: Example of normalized energy distribution functions 
for a = e, t neutrinos at the neutrino sphere in the su- 
pernova environment. Here we take (E Vc ) = 10 MeV and 
(Ev T ) — 27 MeV and the neutrino degeneracy parameter for 
both flavors to be r) Va = 3. 



From Eq. (|78|l it is clear that this term could be sub- 
stantially smaller than Bf^ s if E v is a typical neutrino 
energy. Even if this is not true for E v = E rcs at very 
high density where E res is small, higher energy neutrinos 
and antineutrinos may experience significant in-medium 
mixing angles over a broad range of energy. Though not 
strictly our BDS, this may nevertheless approximate it. 

Previous numerical simulation work on neutrino flavor 
evolution in the supernova environment may offer only 
limited guidance here. The simulation in Ref. [f| made 
the same 2x2, and one-dimensional approximations as 
we make here. (By "one-dimensional" we mean that fla- 
vor histories on neutrino trajectories of any polar angle 
are taken to be the same as a radially directed path for 
the same lapse of Affine parameter along these trajecto- 
ries.) Additionally, the work in Ref. fi} employed the 
density profiles and neutrino fluxes of the Mayle & Wil- 
son late time supernova models and it adopted a range 
for Sm 2 which is now known to be un-physically large for 
active-active neutrino evolution. Both of these features 
combined to produce only minimal effects from rather 
small values of B eT . 

Likewise, the numerical simulation of Ref. consid- 
ered one-dimensional, 2x2 neutrino flavor evolution with 
un-physically large mass-squared difference. The conclu- 
sions in this work regarding real supernovae are suspect 
because: (1) the large Sm 2 used would in reality demand 
the incorporation of sterile neutrinos which mix signif- 



the resonance energy will also increase. It could increase 
significantly if |A + _B| — > 0. If neutrino flavor conversion 
in the channel v e ^ v T is efficient and complete, then 
at some point we will have the situation depicted in Fig. 
Q. Here B could be negative because we have swapped 
flavors at low neutrino energy and, for our chosen spec- 
tral parameters, the v e population now may be smaller 
than the v r population. Furthermore, in this situation 
the material may be driven more neutron-rich (lower Y e ) 
on account of the now altered competition between the 
processes in Eqs. & ©. Eventually, of course, the 
resonance will sweep through the higher energy regions 
of the distribution functions and the fluid element will 
move further out to where neutrino fluxes are lower. 



The higher the resonance energy, the greater the neu- 
trino population which has been appreciably mixed and, 
hence, the larger will be B eT . The importance of this can 
be ascertained by comparing -B^ T DS to the vacuum term 



(86) 
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icantly with actives and this was left out; and (2) the 
feedback of neutrino flavor conversion on Y e was not cor- 
rectly modeled since the threshold was neglected in the 
rate for D e + p — ► n + e + Q and the weak magnetism 
corrections Q were also neglected. However, this nu- 
merical simulation was the first to follow neutrino phases 
in detail in this environment. Synchronization of large 
amplitude neutrino flavor oscillations was seen. This be- 
havior is at least qualitatively like some aspects of the 
BDS, especially as regards significant in-medium mixing. 



Though the conditions for establishment of the BDS 
are manifest in many regions of the post-shock supernova 
environment, it has not been seen unambiguously in sim- 
ulations to date. However, there is considerable room for 
improvement in the sophistication of these simulations. 
Flavor evolution histories on different neutrino trajecto- 
ries needs to be followed in detail, including all coupling. 
The role of density fluctuations [25J in getting some neu- 
trino conversion going at high density also needs to be 
investigated. Likewise, legitimate three-neutrino mixing 
of neutrinos and antineutrinos must be followed. Finally, 
the effects of neutrino mixing on neutrino transport in 
the neutron star core may be important and recent for- 
mulations |26l l2?| of this problem represent significant 
progress. 
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FIG. 4: Neutrino energy distribution functions in a fluid ele- 
ment at some time t corresponding to position r. These have 
the same neutrino energy/temperature parameters as in the 
previous figure. Here, however, complete flavor conversion 
between i/ e and v T has taken place from E v = to the MSW 
resonance energy at this time/position E TCS = 12 MeV. As the 
fluid element moves out the resonance energy will increase and 
sweep from left to right through the neutrino distributions 



C. The Ephemeral Nature of the BDS 

Changing neutrino luminosities and fluxes and chang- 
ing matter density will quickly lead to the development 
of complex amplitudes in the unitary transformation be- 
tween the neutrino mass/energy and flavor bases which, 
in turn, will lead to complex potentials. This will signal 



J 



the end of the strict validity of our particular BDS dis- 
cussed above. However, it may not signal the immediate 
end of appreciable in-medium mixing among the flavors 
of neutrinos and antineutrinos. 

If we ride along with a fluid element being driven from 
the neutron star's surface by heating we will see a lo- 
cal fall off in matter density and neutrino fluxes and 
so a decrease in neutrino-electron and neutrino-neutrino 
forward scattering-induced potentials in this Lagrangian 
frame. What is the effect of this time dependence on in- 
medium flavor mixing? Using the flavor basis evolution 
equation (Eq. I36JI and ignoring the term proportional to 
the identity we can find a second order equation for, e.g., 
a ea , the amplitude for a neutrino of initial flavor a to be 
a v P : 



, 2 _ ^eT . 



(87) 



Here the dots over quantities denote time derivatives and 



\B R 



2iS-2i- 



5B P 



B P 



with S = A + B - A cos 29 and S = A + B. 

In our BDS all time derivatives vanish, u> ~ \B eT \/2, 
and we can solve Eq. (JHZJl, use unitarity (|a eQ ,| 2 + |a rQ | 2 = 
1) and take a ee = a eT = a Te = ± exp (±iu>t)/\/2 and 
a TT = =p exp (Hujt)/y/2, and likewise, a ee — a eT — a Te = 
± exp (zLiojt) I \[2 and a TT = T exp (±iujt)/V2. If we em- 
ploy these solutions in the general flavor-basis form for 
the off-diagonal potential, 



B eT = V2G F / ] (1 — cos6> pq ) [dn Va a ea a* a - dn^ a a ea a* a ], 
we will recover the BDS form for this [c/., Eq. I|77|l ] discussed above: 



Sf T DS w V2G F I (1 - cos6lp q ) [{dn Ve - dn Ve ) - (dn VT - dn^ 



(89) 



(90) 



However, once we allow the potentials to change in time, 
amplitudes will quickly acquire a non-sinusoidal time de- 
pendence which will lead to the development of poten- 
tials with imaginary components. With complex poten- 
tials we will lose a key assumption used in obtaining the 
BDS of Eq. H9U|) . Flavor evolution from that point on 
will be complicated, but there is nothing in the evolution 
equations that demands an immediate return to medium- 
suppressed flavor mixing for most neutrino energies. 



IV. THE BDS IN LEPTON-DEGENERATE 
COSMOLOGIES 



Coherent active-active neutrino flavor evolution in the 
early universe also can be dominated by the flavor off- 
diagonal potential whenever significant net lepton num- 
bers reside in the neutrino seas. Collision- associated de- 
coherence dominates neutrino flavor conversion in the 
early universe at temperatures above Weak Decoupling, 
T > 1 MeV. Neutrino inelastic scattering rates are large 
compared to the expansion rate in that regime. By con- 
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trast, neutrino flavor evolution below this scale is largely 
coherent. We will concentrate on this epoch. For il- 
lustrative purposes we also will confine our discussion 
to cosmologies which have identical lepton numbers and 
interactions for both mu and tau neutrinos. With this 
condition we can reduce the flavor evolution problem to 
the same 2x2 channel we dealt with for the supernova 
environment. 

Coherent medium-enhanced flavor transformation in 
the channel v a ^ i/g (where a, /3 = e, fi, r and a ^ (3) in 
the early universe is governed by the flavor-basis Hamil- 
tonian 



H 



EU 



If y-Acos26» 
2 I A sin 29 + B Te 



A sin 20 + B eT 
A cos 20 - V 



(91) 



where V = A + B + SC, we take a = e and = r, and 
we ignore the components of the coherent Hamiltonian 
proportional to the identity. The time evolution of the 
flavor and mass/energy amplitudes can be handled by a 
mean field Schrodinger-like equation in complete analogy 
to Eq.s & gTJ. 

The thermal contribution to the flavor-diagonal poten- 
tial is SC. The high entropy of the universe (S ~ 10 10 ) 
dictates a sometimes appreciable contribution to neu- 
trino effective mass stemming from forward scattering 
on thermal fluctuations (cf. Ref. |28|). The thermal 
term relative to zero potential is C w r a GpT 5 , where 
r e rj 79.34 and r T « 22.22 at the epoch of interest, 
T < 2 MeV. The difference of these contributions is 
SC w SrGpT 5 , where Sr = r e — r T . This term in 
the Hamiltonian is negligible in the post Weak Decou- 
pling epoch environment whenever the lepton numbers 
are significantly larger than the baryon-to-photon ratio, 
KfJ ^ V- We therefore will neglect SC here, as we 
consider the large lepton number case only. Similarly, 
we also neglect thermal contributions to the flavor off- 
diagonal potential. 

Let us consider the case where £ Vr > £ Vb . We im- 
pose the BDS for this case by: (1) assuming completely 
adiabatic neutrino flavor evolution; and (2) employing 
the maximal mixing angles of Eq.s l|63|) . (|64[) . (|65fl . (I66|) . 
The first of these conditions allows the use of the adia- 
batic forms for B and B eT , Eq.s and i|55|). respec- 
tively. Employing the second assumption in these poten- 
tial forms immediately implies that B as 0. 

We will have \a\ T \ 2 = 0, |ai e | 2 = 1, |ai r | 2 = 1, and 
|aie| 2 = in the adiabatic limit when £ Vt > £ Ue > 0. 
(There is a level crossing for the antineutrinos in this case 
but not for the neutrinos.) Employing these amplitudes 
in Eq. (|55|) and noting the isotropic nature of the neutrino 
distribution functions in the early universe, we find 



B ET « s/2G F [{n Ur -np T )- (n„ e - n p J] (92) 
„ 2^2C(3) GfT 3 



7T^ 



For this solution to be self consistent, we must have 
\B eT \ ^> A and adiabatic neutrino flavor evolution. The 



first condition will be true so long as 

(4 T -4.)»n»7~3xlO- 



10 



(94) 



This condition follows on noting that A = \/2G-pYen\, 
and the baryon density is — r\ n 7 . 

Adiabatic flavor evolution in this case is all but guaran- 
teed on account of the slow expansion rate at this epoch 
and the known values for the atmospheric and solar mass- 
squared difference scales (see the discussion of adiabatic- 
ity in Ref. |29j). The effective density scale height for 
weak charge at an epoch with temperature T in the early 
universe is roughly a third of the causal horizon length 
or, more correctly, 



n = 



1 dV 


-l 










V~dt 




3 



g/g A(/A£ 
~3H 3H 



(95) 



The statistical weight in relativistic particles at this 
epoch is g and the difference in lepton numbers here is 
At = £ Ut — £„ c . The expansion rate (Hubble parameter) 

is H w (87r 3 /90) 1/2 3 1 / 2 T 2 /mpi and the horizon length 
in dff(t) =2t = H^ 1 . Neglecting the rate of change of g 
and the rate of change of lepton difference, At, relative 
to the expansion rate, the adiabaticity parameter for the 
BDS [cf. Eq. (JH21] is in this case 



(96) 
(97) 
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In the second approximation we have used the measured 
baryon-to-photon ratio r\ rs 6 x 10~ 10 . For an epoch 
where T ~ 1 MeV, we will have g « 10.75 and Y e as 
0.5 to 0.16, so that we have adiabatic neutrino flavor 
evolution for A£ ^> r\. Since this is the same condition 
as that required for \B er \ > A to obtain, we see that 
B^ r DS indeed can be a self consistent solution in the early 
universe. 

If we had a lepton number degenerate cosmology with 
£v T > £u e as described above, then we would expect 
a rapid "evening up" of the electron and tau lepton 
numbers if the maximal mixing BDS obtained. In fact, 
numerical simulations which include the neutrino back- 
ground terms and follow neutrino amplitudes with phases 
find just this end result 0. They also find that neutrino 
oscillations are synchronized in phase in time/space. To 
what extent is this result related to the BDS? 

In contrast to the supernova case, it may be much 
easier to evolve into the BDS in a large lepton number 
early universe scenario. First, with modest lepton num- 
bers, \£ Va \ < 10~ 4 , there could be significant transfor- 
mation of neutrino flavors associated with collisions and 
de-coherence in the epoch where temperatures are above 
the weak decoupling scale. The extent of this conversion 
depends on vacuum mixing angle and, hence, possibly 
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on 813. Significant flavor conversion at this epoch would 
imply that there may already be a significant \B eT \ by the 
time neutrino flavor evolution is dominated by coherent 
forward scattering in the epoch after weak decoupling. 

Also, depending on A£, neutrino resonance energies 
could be appreciable at this epoch. The resonance con- 
dition implies that the resonant scaled energy is 



1 



7r 2 Sm 2 cos 26 



4V2C(3)G F (A€ + r?F e ) T 4 



(98) 



It would not take much coherent flavor transformation 
to produce a prodigious B eT , at least on the scale of the 
vacuum term in Eq. I|8tj|) . This term can be especially 
small here on account of the large resonance energies. For 
example, using Eq. (|9"5|l , the resonance energy = Te rcs is 



(6.2 x 10 5 MeV) 



3 x 10" 3 eV 2 
I 



/MeV\ 3 1 

{ T ) Y e + A£/ V 



(99) 



In the BDS only the A term determines resonance ener- 
gies. This case would correspond to A£ — in Eq. iff^)) . 
This would give a large resonance energy and would likely 
lead to a significant \B eT \ and appreciable mixing across 
a broad range of neutrino energies. 

As outlined above, the numerical simulations of Ref.s 
|15| show synchronization of large amplitude neutrino 
and/or antineutrino flavor oscillations in large lepton 
number early universe scenarios. As for the supernova 
case, this looks at least qualitatively similar to the BDS. 
This is especially true while A£ is still large in these 
models. As this quantity is decreased by neutrino fla- 
vor mixing, the flavor off-diagonal potential is decreased 
and neutrino flavor inter-conversion will be less efficient 
and the BDS conditions may no longer apply. 



V. CONCLUSION 

We have investigated the complicated problem of 2 x 2 
coherent neutrino flavor evolution in the limit of large 
flavor off-diagonal neutrino-neutrino forward scattering 
potential both in the early universe and in the post-shock 
supernova environment. We have identified a simple so- 
lution/limit in this problem. This solution (BDS) is gov- 
erned by a dominant off-diagonal potential. This consti- 
tutes a viable solution only under a number of restrictive 
assumptions, but it is evident that even a rough facsimile 
to this solution will retain key principal features of the 
BDS. These include maximal or near maximal in-medium 
mixing angles for both neutrinos and antineutrinos over 
broad ranges of neutrino/ antineutrino energies. These 
features are very different from conventional neutrino fla- 
vor amplitude evolution with the MSW effect. 

Indeed, it has been generally thought that the small 
values of neutrino mass-squared difference among the ac- 
tive neutrinos preclude significant effects from medium- 
enhanced neutrino flavor conversion in the dense environ- 
ment above the neutron star in supernova models. This 
is largely the case for conventional MSW neutrino evo- 



lution. It need not be the case when neutrino-neutrino 
forward scattering potentials are large. Indeed, we have 
outlined above how the neutrino potentials may conspire 
to engineer significant neutrino/antineutrino transforma- 
tion even when Sm 2 values are small. 

We have identified regions and plausible conditions in 
the post-bounce and post-shock supernova environment 
where the BDS, or something approximating it could ob- 
tain. Likewise, we have found that early universe sce- 
narios with significant lepton numbers also provide con- 
ditions favorable for the BDS solution to reign. It is 
not yet clear that either of these venues provides a clear 
and compelling evolutionary path into the BDS regime. 
However, numerical simulations have provided hints that 
something like the BDS may occur in these environments. 
A central question that we leave for the computational 
neutrino flavor evolution community is whether, and/or 
to what extent, the BDS is attained. 

The stakes may be high. If neutrino energy spectra or 
fluxes for the different neutrino flavors are appreciably 
different at any point in the some 20 s time frame follow- 
ing core bounce, then neutrino and antineutrino mixing 
could affect shock re-heating physics, neutrino heating 
feedback on the neutrino-driven wind, and slow outflow 
r-process scenarios. Of course, the neutrino signal could 
be affected by any kind of neutrino/antineutrino flavor 
mixing. The effect we point out here, if it is ever realized 
deep in the supernova envelope, could appreciably al- 
ter the emergent neutrino energy spectra and fluxes over 
those calculated via conventional MSW evolution alone. 

Finally, our considerations extend to any environment 
where neutrino fluxes are appreciable and where neutrino 
flavor mixing may have important consequences for the 
neutron-to-proton ratio and/or energetics and dynam- 
ics. Gamma-ray burst fireball models sited in the vicin- 
ity of a hot or collapsed compact object are a case in 
point. In this environment, the supernova problem, and 
in the early universe we are hard pressed to follow com- 
putationally the evolution of the neutrino/antineutrino 
component. It is unsatisfactory that this remains true 
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even in the face of the tremendous strides in experimen- 
tal neutrino physics which have given us the neutrino 
mass-squared differences and most of the vacuum mixing 
parameters. 
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